Bose condensation in flat bands 
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We derive effective Hamiltonians for lattice bosons with strong geometrical frustration of the 
kinetic energy by projecting the interactions on the flat lowest Bloch band. Specifically, we consider 
the Bose Hubbard model on the one-dimensional sawtooth lattice and the two-dimensional kagome 
lattice. Starting from a strictly local interaction the projection gives rise to effective long-range 
terms stabilizing a supersolid phase at densities above Vc = 1/9 of the kagome lattice. In the 
sawtooth lattice on the other hand we show that the solid order, which exists at the magic filling 
Vc = 1/4, is unstable to further doping. The universal low-energy properties at filling 1/4 + Sv are 
described by the well-known commensurate-incommensurate transition. We support the analytic 
results by detailed numerical calculations using the density-matrix renormalization group and exact 
diagonalization. Finally, we discuss possible realizations of the models using ultracold atoms as well 
as frustrated quantum magnets in high magnetic fields. We compute the momentum distribution and 
the noise correlations, that can be extracted from time of flight experiments or neutron scattering, 
and point to signatures of the unique supersolid phase of the kagome lattice. 

PACS numbers: 75.10.Jm, 67.80.K-, 67.85.-d 
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Strong geometric frustration can prevent straightfor- 
ward ordering and thus lead to the emergence of novel 
highly correlated ground states. The best known exam- 
ples of this phenomenon are from spin systems. Frus- 
tration of the magnetic exchange interactions on certain 
lattices gives rise to extensive degeneracy of classically 
ordered states, ^^^ invalidating a direct semiclassical spin- 
wave analysis. This picture has close analogy in the 
physics of the fractional quantum hall effect, where the 
huge degeneracy of a partially filled Landau level inval- 
idates perturbative analysis in the interactions. In both 
cases the true ground state, which could be a Laughlin 
state, a spin liquid or some unexpected broken symmetry 
state, emerges from the degenerate manifold in a highly 
nontrivial way. 

In this paper we address a related question concern- 
ing the ground states of weakly interacting bosons in 
a lattice which fully frustrates the bosons' kinetic en- 
ergy. The usual expectation is that weakly interacting 
bosons will form a condensate in the lowest-energy single- 
particle state, or in other words, the lowest eigenstate of 
the kinetic-energy operator. However, if the hopping ma- 
trix elements on the lattice are sufficiently frustrated, the 
lowest Bloch band becomes flat, thus providing a huge de- 
generacy of single-particle states to which the bosons may 
condense. The nature of the ground state is now fully 
determined by the interactions acting within the hugely 
degenerate manifold. Under these conditions a straight- 
forward perturbative treatment in the interaction is of 
no use. The problem is inherently strongly correlated 
and provides an interesting route for understanding and 
perhaps even realizing novel phases of matter. 

We shall specifically consider a Hamiltonian of the 



H = Y.\'.Mb^ + M + ^T.^\b\hb. 



(1) 
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where bi are bosonic operators defined on the sites of the 
two-dimensionalkagome lattice. This model gives a flat 
lower Bloch band in the single particle spectrum. ^°"^^ 
We shall also consider a related one-dimensionalmodel 
defined on the sawtooth lattice. Both models and the 
band structure they give rise to are depicted in Fig. 1. 

Such models of bosons with flat bands are of direct 
relevance to real physical systems. Recently a number 
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FIG. 1. (color online) (a) Single particle dispersion on the 
kagome lattice along high symmetry lines in the Brillouin zone 
(gray), (b) The kagome geometry with its lattice vectors ai/2 
and the basis sites A, B, and C in each unit cell (gray), (c) 
Single particle dispersion on the sawtooth lattice as a function 
of momentum k for t' — \/2t. (d) The sawtooth geometry 
with couplings t and t' and the basis sites A and B in the 
unit cell (gray). 
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FIG. 2. (color online) Strictly local eigenstates of the ki- 
netic energy (local magnon states), (a) Resonating hexagon 
state on the kagome lattice where the + and — refer to the 
wave-function amplitude on each site around the hexagon, (b) 
Localized V-state on the sawtooth chain, (c) Ground state of 
the many-magnon system on the kagome lattice at density 
Vc = 1/9. (d) Ground state of the many-magnon system on 
the sawtooth chain at i^c = 1/4. 



of proposals were put forward for realization of models 
with frustrated hopping using ultracold atoms in opti- 
cal lattices. ^■^'^^ Another natural realization involves frus- 
trated spin-1 magnets. If the Curie- Weiss temperature is 
sufficiently low, as in to-MPYNN-BF4 (Refs. 15 and 16) 
(Bcw~3K), the spins can be fully polarized, or nearly 
so, by external magnetic fields. The dilute population 
of magnons, or depolarized spins, in the highly polarized 
regime is well described by Hamiltonian (1). 

The presence of a flat band implies the existence of 
localized eigenstates of the kinetic energy, as illustrated 
in Fig. 2(a) and 2(b) for the kagome and sawtooth lat- 
tices. At sufficiently low filling of the lattice one can 
construct exact many-body ground states of Hamiltonian 
(1) from the localized single particle states. ^^^^^'^'^ This is 
done by occupying some subset of spatially nonoverlap- 
ping localized states. Of course the ground state in this 
case is massively degenerate, because of the many pos- 
sible arrangements of nonoverlapping localized orbitals. 
However, there is a critical density Vc at which the oc- 
cupied localized states are closely packed, making an or- 
dered crystal structure that breaks the lattice symmetry 
as shown in Fig. 2(c,d). The critical density is Vc = 1/9 
for the kagome and Vc = 1/4 for the sawtooth lattice, re- 
spectively. Formation of the close packed crystals is ex- 
pected to give rise to magnetization plateaux in quantum 
magnets^^ and to incompressible insulators in ultracold 
atomic systems. Any additional bosons on top of the 
close packed configurations must overlap spatially with 
other bosons and the exact construction then no longer 
works. 

The question we address in this paper is what ground 
states form at lattice filling slightly above Vc- In partic- 
ular, do the added atoms condense to form a superfluid 
on top of the underlying density wave? Such condensa- 
tion is unusual because it would be entirely driven by the 
interactions rather than by the hopping. Another ques- 
tion is whether delocalization of the added atoms leads 
to destruction of the solid order. 

Our general strategy to treat these systems is to 



project Hamiltonian (1) onto the flat Bloch band and to 
obtain in this way an effective low-energy Hamiltonian 
which depends only on the (weak) interaction U . This 
approach is in the same spirit in which projection to the 
lowest Landau level is used to derive effective theories of 
quantum Hall states. ^"^ The resulting low-energy model 
is defined on a new lattice, and is not frustrated. We are 
therefore able to analyze it using standard methods, such 
as bosonization and mean-field theories. The results are 
then checked against numerical calculations performed 
on the original Hamiltonian (1) defined on the kagome 
and sawtooth lattices. 

Before proceeding we make note on previous work 
dealing with lattice boson models such as Eq. (1) in 
the regime of large number of bosons per site i^ 3> 1, 
which is the opposite of the regime we consider here.^^ If 
V ^ \t\/U > 1, then Hamiltonian (1) can be mapped to 
an effective Josephson array model. In the classical limit 
(where the charging energy can be neglected) the ground 
state manifold can be mapped onto a height model at 
its roughening transition. As a consequence, there is no 
long-range order in the phase (e*"^), but a "tripled" con- 
densate forms where {e^'^'^) exhibits long-range order. ^-'^ 

Below we lay out the structure of the paper and give 
a brief overview of our main results. In Sec. II we con- 
sider the sawtooth model. The one-dimensional char- 
acter of the chain allows to introduce the key ideas in 
a simple setup before applying them to the more com- 
plicated kagome lattice. The crucial step is to define 
an orthonormal basis of Wannier functions, which reside 
in the fiat band. Projection of the Hamiltonian on the 
fiat band then takes the form of a spin-1/2 model with 
the spins located on the center of the Wannier functions 
and representing occupation of these states by bosons. 
Next, we take the long wave-length limit in a harmonic 
fiuid approach. ^^ This analysis shows that upon increas- 
ing the density away from the magic filling i^c = 1/4, 
the charge density wave is immediately destroyed due 
to proliferation of domain walls, which gain their mobil- 
ity from the interaction term. The transition is in the 
universality class of the commensurate-incommensurate 
transition, ^"^"^^ from which several sharp predictions fol- 
low. First, the delta- function peak in the density corre- 
lations at wave-vector tt/o, which appears at 1/4 filling, 
splits into a pair of power-law singularities at ('K±2Si>)/a, 
where Si^ = v — 1/4 > 0. Second, the established 
condensate is characterized by the Luttinger paranreter 
K — {1 + I?(5i/)/4, which controls the power-law singu- 
larities. Here I? is a positive constant independent of v. 
Note, that the value oi K = 1/4 at the transition point 
is quite remarkable for a system of bosons with purely 
on-site interactions and is a direct consequence of the ge- 
ometric frustration. We find all of these predictions to be 
in good agreement with numerical density-matrix renor- 
malization group (DMRG) calculations performed on the 
original sawtooth model. 

In Sec. HI we turn to the kagome lattice, focusing on 
filling ly > Vc = 1/9. Added bosons may either form mo- 



bile domain walls in the solid order or alternatively move 
as independent particles on interstitial sites of the den- 
sity wave. Within an effective theory for these defects we 
estimate that the energy of a mobile interstitial, which 
does not destroy the density wave order, is lower. To 
substantiate this result we again derive an effective low- 
energy model by projecting Hamiltonian (1) onto the flat 
band. The resulting spin-1/2 model, written in terms of 
occupation of the flat band Wannier states, lives on the 
triangular lattice defined by the mid points of hexagons 
in the original kagome lattice. Because of the absence 
of a gap in the single particle spectrum [see Fig. 1(a)], 
the justification of the low-energy theory is not as obvi- 
ous in this case as in the sawtooth model. However, we 
show that interaction with the uniform component of the 
density in the flat band generates a gap to excitation of 
the second band, which justifies the projection scheme. 
A mean-field treatment of the low-energy Hamiltonian 
yields a "lattice supersolid" phase in the intermediate 
density regime Vc < v < Vm-e, in which both the charge 
density wave (CDW) order and the off-diagonal order are 
present. The phase ordering structure is in the form of 
a 120° pattern on the a/3x V3 unit-cell of the CDW. In 
Sec. IV we discuss the experimental implications of our 
finding for both the cold-atoms and the spin-1 implemen- 
tations. 



II. SAWTOOTH CHAIN 

We start with the discussion of the sawtooth chain 
model^'''^^'^^ [cf. Fig. 1(d)], which can be viewed as a 
one-dimensional analogue of the kagome lattice. It is 
similarly constructed from corner-sharing triangles and 
can support a flat band for a specific ratio between the 
hopping matrix elements along the base and edges of the 
triangles. 



A. Flat band 

Let us review the features of a tight-binding model 
on the sawtooth lattice. We shall obtain the ratio of 
coupling constants required to give a flat lower band and 
define basis sets of localized states residing in this band. 
This will be important for implementing the projection 
on the flat band in the next section. 

The hopping Hamiltonian can be formulated as a 2 x 2 



matrix in terms of the "spinor" 6^ — [fos./t, &A,fc] 
senting the two sites in the unit cell: 



H,, 



Eh 



2tcos{ka) t'il + e''"'') 



repre- 



(2) 



The sum runs over the first Brillouin zone and a denotes 
the distance between adjacent B-sites. The resulting dis- 



persion is given by 

huj±{k) = tcos(fca) ± ^/f^ cos2(fca) + 2i'2 cos(A:a) + 2t'^ 

t''W2t J — 2t, 

~ I 2t[l + cos(fca)]. 



(3) 



We see that for the special ratio 7 = t' /t = \J'l we indeed 
obtain a flat lower band which is separated by an energy 
gap A = 2t from the next band. The eigenstates take 
the form 






-sin(t?fe/2) cos(z?fc/2)e-*'="/2 
cos(i?fe/2) sin(i5fe/2)e^'='*/2 



(4) 



i?fc = arctan 



2^/1 + cos(fca) 
cos(fca) 



7r8[— cos(fca)]. 



The flat dispersion allows to construct wave-packets of 
Bloch states (4) from the flat band that form localized 
kinetic energy eigenstates. For example, we can obtain 
the states shown in Fig. 2 that are strictly localized to 
three sites, as 



Vi 



vjio; 



^(^/26l,, 






Kd\^) (5) 



Because these states are strictly localized, one can use 
them to construct many-body ground states of the in- 
teracting Hamiltonian by occupying only nonovcrlapping 
] Vi ) states. In particular the CDW state on a lattice 
with N sites is given by 



Af/2 

|CDWi/4)-nvL 



(6) 
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This state is clearly an eigenstate of the interaction en- 
ergy with zero eigenvalue as there are no two sites with 
more than one particle. Note that A^/2 particles corre- 
spond to a density v^ — 1/4. 

Although the states ] Vi ) are linearly independent and 
complete in the flat band subspace, they are not a conve- 
nient basis to work with because they are not orthogonal. 
This is evident from the fact that neighboring states ] Vi ) 
and ] Vij^x ) actually share a site. We therefore define an 
alternative set of localized wave functions residing in the 
fiat band, which do form an orthonormal basis. These 
are constructed as Wannier states. ^^ 

W^/ I 0) = E {""A^^^ - ^^)bh + '^*Birj - r.)6L,j) I 0) , 



WB{rj) 



U2^ 



"" dk 



sm{^k/2)e'''^ 



WAirj)= I |^cos(^,/2)e''='^^"^^/^ 



(7) 



(8) 



By construction these states are orthogonal and expo- 
nentially localized around the unit cell i. Locally, they 
have the structure of the V-states but they fall off expo- 
nentially with a localization length of ^ ~ log(2.15)a, cf. 
App. A. 



B. Projection onto the fiat band 

The exact ground state ICDW1/4) at Uc — 1/4 cor- 
responds to a close-packing of V-states. If we dope the 
system with additional particles we loose this property 
and we expect a transition to a different phase. Let us 
assume small interactions [/ < A. Any phase transition 
as a function of v is then necessarily driven by the in- 
teractions as the kinetic energy is completely quenched 
for energies below A. This provides strong motivation 
for projecting the interaction Hamiltonian onto the flat 
band, in analogy with projection to the lowest Landau- 
level for quantum Hall states.^" 

To implement the projection we express the bosonic 
operators in terms of the Wannier basis derived above 



^A{B),i - Z^^A(B) 



(r,j — rj)W^ + higher bands. (9) 



iour-W operator terms. Even off-diagonal terms exclu- 
sively include assisted hopping ex /■*"(=?). The processes 
in Hamiltonian (10) are illustrated in Fig. 3. 

Note that Hamiltonian (10) describes a spin system on 
a linear chain. The sawtooth geometry is absorbed by 
the projection onto the flat band. The density i^c of the 
original model now translates to half-filling, or zero mag- 
netization of the effective spin-1/2 model. Furthermore, 
the exact CDW solution of the underlying model corre- 
sponds to a perfect Neel order of the spin chain. Due 
to the truncation of the Wannier wave functions, this 
state is not an exact eigenstate anymore. However, the 
dominating Ising interaction I^ 3> /''"("p' still leads to a 
large staggered magnetization in the ground state at zero 
magnetization. 



C. Analysis of the effective Hamiltonian 



And project out the contributions from the higher band. 
Due to the exponential tail of WA{B)i'r) the projected 
Hamiltonian acquires arbitrarily long-ranged interaction 
terms. However, the rapid decay of the Wannier func- 
tions on the scale of a lattice constant allows us to trun- 
cate the sum to \ri — rj\ < a. 

The resulting Hamiltonian contains a renormalized on- 
site interaction (U) and further range interactions as well 
as assisted hopping terms. Since all the coupling con- 
stants turn out to be smaller than the on site repulsion 
by a factor of ^ 4 we replace U with a hard core con- 
straint. Now the effective Hamiltonian can be written as 
a spin-1/2 chain model: 

i/.„ = J2 I'<^t<^Ui + ^ i^t + 1/2) (a+±iar^2 + H.c) 

i 

+ /- [{at + 1/2) (a+ lar^i + H.c)] . (10) 

The interaction constants derive from the overlap of 
Wannier functions on different sites and are given by (see 
App. A for more details) 



We are now in the position to investigate the fate 
of the CDW at lattice filling above Vc, which corre- 
sponds to (cr^) > in the effective model. The Hamil- 
tonian (10) is similar to a AAZ-model with a dominat- 
ing Ising interaction, but with hopping terms replaced 
by assisted hopping. In the latter model the charge 
density is destroyed upon doping away from the magic 
filling by the mechanism known as the commensurate- 
incommensurate transition. ^■^"^''^ We will show that the 
system described by Hamiltonian (10) undergoes the 
same transition. 

As usual in one dimension single particles doped into 
the commensurate CDW can either go into an interstitial 
site of the CDW or break up into a pair of domain-walls 
in that order, each carrying half of the added charge.^® 
We will show below, that energetically, condensation of 
domain walls is favorable to that of the interstitial sites. 
With this input we shall derive the appropriate long 
wavelength description of (10) via bosonization, which 
yields specific predictions for the emergent critical state 
at v > Vc- 



U^ 0.40C/(^oo), 
I"- w -0.025 [/, 



r« 0.112 L/, (11) 

r^ « -0.011 u. (12) 



Note that all terms are proportional to the interac- 
tion U. This is because the effective Hamiltonian (10) is 
just the projection of the interaction onto the fiat band. 
For the same reason all terms in the Hamiltonian are 
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FIG. 3. Processes in the spin-1/2 Hamiltonian (10). Arrows 
denote hopping (or spin-exchange) and dashed lines stand for 
density (a^) interactions. 



1. Domain walls versus interstitial condensate 

Particles added to the CDW state correspond to 
adding up-spins to the effective spin model on top of the 
z-Neel state. The flipped spins can behave in two ways: 
(i) they can exist as well defined excitations and condense 
as interstitials without destroying the CDW order; or (ii) 
break up into two mobile domain walls in the Ising spin 
configuration, each carrying a spin-1/2, thus leading to 
destruction of the CDW order. If one thinks of domain 
walls as the elementary excitations, then the option (i) 
above can be viewed as having a bound state of a kink 
and an anti-kink. Which of the two options is realized in 
a specific model is ultimately a matter of energetics. 

Let us then estimate and compare the energies of an 
interstitial to that of a pair of free domain walls. The 



Ising interaction energy associated with an interstitial 
(flipped spin) is 2/^ because this state violates two Ising 
bonds. In addition, the interstitial can hop to the next 
empty site (down-spin) via the term P^. Therefore, the 
kinetic energy gain is — 2/"^. 

The pair of domain walls still violate two bonds and 
therefore incur the same interaction energy as the in- 
terstitial, i.e., 2/^. However, each domain wall can move 
due to the action of the assisted hopping term J"^", which 
leads to the kinetic energy gain —4/'^" per added particle 
(two domain walls). With the coupling constants given 
by (12), the pair of domain walls are seen to be clearly 
favorable to the interstitial. We therefore conclude that 
the long-range order of the CDW state at 1/4 filling will 
be immediately destroyed by proliferation of free domain 
walls upon increasing the density. 



2. Long wavelength description 

The energetic considerations detailed above imply de- 
struction of the density wave state immediately upon 
increasing the density above v^ = 1/4. The emergent 
critical state is naturally described within the contin- 
uum limit by bosonizing the effective spin-chain model 
(10). However, the usual scheme, ^^ of first mapping to 
a fermion system using a Jordan- Wigner transformation 
is not straightforward. Because there are no quadratic 
terms in the spin-1/2 Hamiltonian, there results no obvi- 
ous Fermi surface to support an expansion in slow modes. 
Under such conditions it is more convenient to employ the 
so called phenomenological approach to bosonization of 
bosons. '^'^ For this, we reintroduce the on-site interaction 
U and write the density and field operators as 
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^i[TT'ni—2'mi} [xi)] 



Note, that we expand around v = l/2a. The most rele- 
vant density wave captured with this approach is of wave- 
length \ — 2a and is described by the order parameter 
(cos(4i?(a:))). We take a naive continuum limit by rewrit- 
ing (10) in terms of the new fields '^{x) and i?(a;), keeping 
only the most relevant terms we obtain 
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(13) 
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We find the expected CDW reproduced also within this 
long-wavelength theory, as the term ex cos(4??) is relevant 
for K < 1/2. Note, that instead of taking the naive con- 
tinuum limit one can regard (13) as a phenomenological 
description of the CDW, as long as K < 1/2. The results 
we draw from (13) will not depend on these details but 
rely solely on the form of the cosine term. 

Let us now proceed with the discussion of an additional 
density Sv — i^ — i'c above the CDW. Forcing Sv > is 
equivalent to a chemical potential larger then the CDW 
gap AcDw To characterize the resulting gapless state 
above Acdw, we can use results from the commensurate- 
incommensurate transition introduced by Pokrovsky and 
Talapov,^^ in particular the critical behavior discussed 
by Schulz.^^ The long- wavelength theory is now given by 
a Luttinger liquid with^^ 



K^^{l+V6iy), 



(14) 



where 2? is a positive number independent of Sv. Further- 
more, we can calculate correlation functions by using the 
new Luttinger liquid with K and by introducing a shift 
'0{x) — ?> dlx) + TTxSiy. One readily obtains^^'^^ 



(p(x)p(0)),«^/,«cos(7r(l + 2Ha;/a)(-J , (15) 
{WHx)Wm,^o^{^y^"' ■ (16) 



We show a sketch of the density correlations in Fig. 4(a), 
together with the v = v^ result. 

From the above considerations we conclude that under 
particle doping the exact ground state ICDW1/4) melts 
due to delocalization of domain walls. The harmonic fluid 
approach allows for a prediction of both the functional 
form of the density-density correlations and the effective 
Luttinger liquid parameter K . In the next section we test 
these predictions by numerical analysis of the sawtooth 
chain. 



D. Numerical analysis 

To complement the analytical study and confirm its 
predictions we analyze the original sawtooth lattice 
model numerically using the DMRG method. Calcula- 
tions are performed on an 80 site lattice [L = 40 unit 
cells) with the parameters t = U = 1 and t' — \/2 and we 
retain up to 200 states. In all calculations we make use of 
the ALPS DMRG code.^^ We consider particle numbers 
iV = 20 -f (5iV ((5iV = 0,1,2,3,4) particles. Note, that 
(5iV = corresponds \,o v = v^^ and we indeed find the 
exact ground-state ICDW1/4) at this point. This is evi- 
dent from the perfect staggered density correlations seen 
in Fig. 4(b) as a sharp peak in the static structure factor 
at wave- vector tt/o. As we dope the system above v^ the 
peak splits in two, with a shift that grows by 27r/L per 
added particle. Exactly as expected from the bosoniza- 
tion study above. 
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FIG. 4. (color online). Fourier transform of the density cor- 
relations in the sawtooth chain, (a) Prediction of the long 
wavelength theory. The delta function peak at g = Tr/a (red 
dashed line), which signals the CDW at i/ = i-c, splits with 
increased density ly = i/c + Sv into two power-law peaks (blue 
full line) at q± = 7r(l ± 2Sv). The effective Luttinger param- 
eter, which controls the power-law is given in Eq. (14). (b) 
Density correlations computed using DMRG. The splitting of 
the peak by one quanta Sk — ^ per added particle confirms 
the predicted q± of the long-wavelength theory. 



To capture the critical behavior of the system it is bet- 
ter to extract the ofF-diagonal correlations, because they 
have a faster decay. At the special point v = i^^ these are 
found to be strictly local, as expected. For SN > on the 
other hand the decay of the off-diagonal correlations is 
consistent with a power-law. As the density is increased 
the critical exponent is reduced continuously from a~ 2 
at 5N — 0+. This is consistent with a Luttinger pa- 
rameter K which grows continuously from K « 1/4, as 
predicted by the long-wavelength theory above, cf. Fig. 5. 

The numerical study thus confirms our main predic- 
tions. The density wave state is immediately destroyed 
upon increase of the density. The exact density wave 
state is replaced by power-law density correlations at the 
incommensurate wave-vector consistent with the doping 
away from Vc- The emergent Luttinger liquid phase that 
forms is characterized by a Luttinger parameter K > 1/4. 



III. KAGOME LATTICE 

We now turn to the discussion of the kagome lattice, 
focusing on the question, which ground state forms when 
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FIG. 5. (color online) The effective Luttinger parameter K 
as a function of density Si/ = i/ — Uc extracted from the off- 
diagonal correlation function calculated using DMRG. The 
black line is a guide to the eye. 



the density is increased beyond i^c ~ 1/9. 



A. Flat band 

Let us review the essential features of the tight-binding 
Hamiltonian on this lattice. It is useful to introduce the 
vector notation 




and 6k 




^e»''-^-6,. 



The unit cells i are at positions Ri = rrnai + nia.2, 
niijUi S Z with ai = a(l,0), 3.2 — a(l/2, -\/3/2), and 
a the lattice constant. Furthermore, the subscripts A, 
B, and C label the lattice sites in a unit cell at Rj, 
Ri -I- ai/2, and R^ -I- a2/2, respectively, cf. Fig. 1(b); 
M denotes the total number of unit cells. The wave vec- 
tors are confined to the first Brillouin zone as defined by 
the Wigner-Seitz cell with respect to the reciprocal lat- 
tice vectors Qi = ^(l,-l/%/3) and Qa = ^(0,2/V^), 
see insets of Fig. 1(a). The hopping Hamiltonian now 
reads 
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where fc,^ = k • a^, j/ = 1, 2 and ^3 ~ ki — k2- Diagonal- 
izing this matrix gives the three bands: 



Swo(k) = -2t, ;ia;±(k) =t[l± v/3T2A(k)], with 
A(k) = cos(fci) + cos(fc3) + cos(fc2), (18) 

shown along the high symmetry lines in Fig. 1(a). The 
eigenstates corresponding to the flat band are given by 
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sin(fci) . 



The existence of a flat band allows to construct local- 
ized eigenstates of the same energy. The most local wave- 
function of this type is confined within a single hexagon^^ 
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Here ai labels the sites around the hexagon, cf. Fig. 2(a) 
and i labels the hexagon. It is easy to see that destruc- 
tive interference on the sites adjacent to the hexagon 
sites prevent the spread of the wave-function. Delocal- 
ization on the hexagon gives rise to the kinetic energy 
—2t. Therefore Al indeed creates an eigenstate with the 
energy of the flat band. Moreover, by acting on the vac- 
uum with the operators A^^ to create any number of spa- 
tially nonoverlapping states, it is possible to the construct 
exact many-body ground-states of the total interacting 
Hamiltonian (1). For sufficiently low lattice filling there 
is a huge number of degenerate arrangements of the lo- 
calized states. However, these localized orbitals become 
close packed at the critical filling ly^. — 1/9, at which point 
there are only three ground states corresponding to the 
three degenerate density wave states: 
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The set T contains close packed nonoverlapping hexagons 
as depicted in Fig. 2(c). 

At higher lattice filling ly > 1/9 such a construction 
would fail to create eigenstates of the interaction term 
because it is no longer possible to construct a spatially 
nonoverlapping set of states. Our strategy to address this 
regime is, the same as in the sawtooth model, to derive an 
effective Hamiltonian using a projection on the fiat band. 
For this purpose we shall introduce an orthonormal set 
of localized states, which belong exclusively to the flat 
band, and are defined on every unit cell. These are the 
Wannier states 



W,{r) 
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where the integral runs over the first Brillouin zone and 
vq denotes its volume. Note, that due to the band- 
touching at the F-point, these Wannier states are not 
exponentially localized but have power-law tails. ^^ Their 
slowest decay is along high-symmetry directions of the 
lattice where they fall off W{r) ~ l/|r|; see App. A for 
more details. 



B. Projection into the flat band 

We are now in the position to derive an effective Hamil- 
tonian for the low-energy dynamics by projecting the in- 



teractions onto the flat band. In this case the projection 
is fraught with further complication because of the van- 
ishing gap in the single particle spectrum between the 
flat band and the second Bloch band. Recall that in 
the case of the sawtooth model the projection was con- 
trolled by the band gap A = 2t, which gave rise to the 
small parameter U/A. The situation is different on the 
kagome lattice because the non-interacting model has a 
topologically protected^^''^^ band touching (A = 0) at the 
F point. Below we argue that for filling factors i^ ^ 1/9 
the interactions open an effective band gap which ulti- 
mately controls the projection. ^^ 

The key point is that a particle in the second Bloch 
band, with small quasi-momcntum k ~ 0, suffers an en- 
ergy shift A w C//9 due to interaction with the mean 
density (i^)j. A particle on the flat band on the other 
hand, does not suffer the same mean-field energy shift, 
because of the extra freedom it has to be in superposi- 
tions of Bloch states. This is explicitly evident in the en- 
ergy scales of the effective spin Hamiltonian generated by 
projection to the fiat band, which turn out to be ^ 0.02C/, 
much less than the energy shift A « U/9 of the second 
band. Thus the projection is justified a posteriori. 

Technically we carry out the projection by expressing 
the boson operators appearing in the original Hamilto- 
nian (1) using the Wannier operators Wi while truncat- 
ing the contribution from Wannier states of the higher 
bands. The resulting Hamiltonian then contains terms of 
the form lijrsWjW^WrWs, where the coefficients lijrs 
consist of wave-function overlaps of the Wannier func- 
tions located at i,j,r, and s. The largest term is an 
effective on-site interaction U « 0.14 [/, which comes 
from i ^ j = r = s. The next coefficient stems 
from i — j and r = s with i and r nearest neighbors. 
The resulting nearest-neighbor interaction has magni- 
tude Iz « 0.02 U <^ U. Other terms, as we shall see 
below are even smaller. We can therefore safely replace 
U with a hard core constraint. This restriction is also re- 
quired in order to be consistent with neglecting occupa- 
tion of particles in the second band, which is at an energy 
A « U/9 above the low-energy manifold, but comparable 
to the on-site interaction. The projected Hamiltonian is 
therefore given by a hard-core boson (or spin-1/2) model 
on a triangular lattice. 

The effective Hamiltonian contains hopping and inter- 
action terms which fall off with range as l/|rp, or faster, 
due to the algebraic decay of the Wannier functions. To 
keep the model simple, we truncate all interaction terms 
smaller than /^„t = 7^/20. When necessary, we check 
against convergence with respect to this cut-off. 

To explicitly write the effective Hamiltonian we de- 
fine the vectors a.a, a — 1,..,6, connecting the nearest 
neighbors on the triangular lattice as shown in Fig. 6. 
The Hamiltonian in terms of spin-1/2 operators CTj , af 
is given by 
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The coefficients read 
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(24) 



In Fig. 6 we illustrate the processes in iJ„ff. 

Note that by projecting to a single band we removed 
2/3 of the degrees of freedom and thus changed the lattice 
structure to a triangular lattice with one instead of three 
sites per unit cell. The critical density at which the den- 
sity wave states form in the kagome lattice is translated 
accordingly to 1/3 filling of the effective model on the 
triangular lattice. Because the repulsive (or Ising anti- 
ferromagnetic) interaction I^ is the largest energy scale in 
(23) we do indeed expect the density-wave state to form 
at this filling. Note, however that the state with exactly 
one boson on every three sites is not an exact eigenstate 
of the effective model (23). This is of course due to the 
truncation of longer range interactions in (23). In what 
follows we will study the properties of the effective model 
and the fate of the CDW order for filling v > 1/3 of the 
triangular lattice. 



C. Domain walls versus interstitial condensate 

When we increase the lattice filling beyond Vc = 1/9 
of the kagome lattice or 1/3 of the effective triangular 
lattice we necessarily introduce defects to the density 
wave structure with one particle in every three unit cells. 




FIG. 6. (color online) Illustration of the processes appearing 
in the effective Hamiltonian (23) for the kagome lattice. A 
dot represents a density operator, hnes denote either a near- 
est neighbor interaction [P) or a hopping (J^'') assisted by 
the nearby density operator. Note that /J'/g have negative 
effective hopping amphtudes (red) while Jj/^ have positive 
hopping (blue). 



There are different types of defects that can accommo- 
date the extra particles. On the one hand, the extra 
particles can hop around interstitial sites of the lattice, 
which is innocuous to the underlying density wave struc- 
ture. On the other hand, they can nucleate mobile do- 
main walls in the CDW order, that ultimately destroy 
it. 

An extra particle in an interstitial site has three nearest 
neighbors in the CDW background. The resulting inter- 
action energy per particle is 3/^. The particle can also 
hop between interstitial sites and thereby gain kinetic en- 
ergy from the assisted hopping matrix elements I^^ and 
I^^ . The delocalization energy of the added particle is 



3|/ri 



-2.5/, 



-6|/f^| 

Let us compare to the case when particles are added as 
domain walls. Domain walls between the three degener- 
ate CDW states come in two flavors, the "straight" and 
"zigzag" domain walls, which are illustrated in Fig. 7. 
The line density, or number of added particles, associ- 
ated with the "zigzag" domain wall is twice that of the 
"straight" domain wall. However, they both incur the 
interaction energy cost of 3/^ per added particle, which 
is the same as the interaction energy of a particle in an 
interstitial site. 

The differences in energy between the defects stems 
from the delocalization energy of particles within them. 
Consider first the "zigzag" domain wall. Half of the par- 
ticles in it can resonate independently on the horizontal 
bonds, as shown by the ellipses in Fig. 7, facilitated by 
the assisted hopping process l2^ . The energy gain per 
added particle is then seen to be — 6|/2^| ~ I^- Particles 
in the "straight" domain wall, on the other hand, can 
hop with the matrix-element /j^ only to reach a point 
where the interaction energy is increased by one bond 
Iz- Therefore, they can only move coherently by a sec- 
ond order process, gaining them delocalization energy of 
order {12^)"^ /Iz- In both cases, of "zigzag" and "straight" 
domain walls, the energy gain from delocalization is less 
than that of the interstitial defect. 

We conclude that doping of particles into interstitial 
sites of the CDW, due to the high mobility they ac- 
quire, is energetically favorable to nucleating domain 
walls. Since the interstitial particles do not destroy the 
CDW, we expect that order to survive and co-exist with 
a condensate formed by the mobile interstitial particles. 
In the next section we use the effective spin model on the 
triangular lattice to formulate a mean-field description of 
the resulting supersolid phase. 




CDW + SF 



FIG. 7. (color online) Three patches of the CDW (different 
colors/gray scale) with two types of domain walls: "straight" 
(left) and "zigzag" (right). The thick lines denote the vio- 
lated bonds which lead to the interaction energy cost. The 
yellow (shaded) ellipses around half of the particles in the 
"zigzag" domain wall indicate resonating particles which gain 
the delocalization /j^. 



D. Mean-field treatment 
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FIG. 8. (color online) Order parameters in the kagome sys- 
tem as a function of the density i/ — Vc + Sf. The lines are 
the result of the mean-field analysis of the effective Hamilto- 
nian on the flat band. The CDW order parameter is shown 
as a full (red) line while the superfluid order parameter is the 
dashed (blue) line. Circles (CDW) and triangles (superfluid) 
mark the values obtained from exact diagonalization of the 
original kagome system on a small cluster. The (blue) star at 
Sf = 0.15 is the fluctuation reduced superfluid order parame- 
ter. The mean-field result shows a critical density where the 
CDW order vanishes. 



In the last section we established that the added par- 
ticles preferably hop on interstitial sites of the CDW 
and do not lead to proliferation of domain walls. Con- 
sequently the CDW order is preserved and a mean-field 
treatment of the coexisting CDW and interstitial conden- 
sates is justified. The mean-field approximation consists 
of making a site factorizable ansatz for the wave-function: 



I*) = [][cos(z?,/2)| i), + e'^^ sin(7?,/2)| t)^ 



(25) 



In order to account for the broken translational symme- 
try, we assume space dependent coefficients "di and tpi 
with the unit cell of the CDW at filling v^ = 1/9. The 
variational energy _E^„ = {'i'\H^ff\'^), then depends on 
the three angles ??q and phases tpa as 
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The sum runs over cyclic permutations of A, B, and C 
and ipAB = ^A — fB, etc. Here we only gave the explicit 
expression for the expectation values of the leading terms 
in (23). When solving the variational problem numeri- 
cally, we also included, and found no significant effect, of 
the other terms. 

Let us first consider two simple limits of the variational 
problem. The pure CDW state corresponds to the con- 
figuration -ffji = TT and ^B = ^c =0. Plugging this into 
(26) we see that all terms, "interaction" and "kinetic". 



vanish and the phases ipa all drop out of the problem. 
The other simple limit is a pure condensate with no den- 
sity modulation, so that -(9^ = i?^ = t^c = i^- In this 
case, the problem reduces to 



E[i3,ipa] =2sin2-<{ {-3\I^y\smi3 + 3h{l - cos{i))) 



+ \I^y\smd[cos{ipAB) + cos{(pBc) + cos(^ca)] 



The value of d is set by the density 5i^ and we are left 
with an optimization with respect to the phases (pa ■ The 
energy is minimal for 

(PA = 0, ^B = ±27r/3, </7c = ±47r/3. (27) 

For the experimental realizations discussed below, it is 
interesting to note that when Fourier transfornred, this 
phase pattern corresponds to the K or K' points of the 
Brillouin zone of the triangular lattice. However, in the 
reduced CDW Brillouin zone these points correspond to 
the r point. 

In the general case we optimize for the angles da and 
phases (pa numerically to find the mean-field density pat- 
tern. Starting from an initial guess of the form ^^ = tt, 
"Sb/c = Oi the minimal energy is always obtained for a 
solution of the form -Oa < tt, ^b = "dc ^ 0. Furthermore, 
above the critical filling (i.e. Siy > 0) the phases quickly 
approach the 120° pattern found above for the case of a 
pure condensate. 

Using the optimized variational wave-function we com- 
pute the order parameters of the CDW 
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and of the superfluid 



^ = matin- 



(29) 



Both are shown in Fig. 8. The CDW order parameter 
I'stag is seen to decrease continuously from the perfect 
CDW at the magic fiUing i^c = 1/9 (Sf = 0) until it 
vanishes at an upper critical density Svj^p sa 0.13. At 
the same time there is a build up of the superfluid order 
parameter ip. The change from an increase to a reduc- 
tion and eventual vanishing of ip with density at t/ = 1/3 
is an artifact of the mapping to spin-1/2 degrees of free- 
dom (hard core bosons) . However at this stage our lowest 
band approximation breaks down. The energy cost of in- 
teraction with neighboring particles grows to a point that 
it becomes energetically favorable to occupy the second 
band (and consequently also induce double occupancy). 

The above analysis of the effective model should give 
qualitatively correct results at moderate densities, but 
should not be relied upon for high quantitative preci- 
sion. One cause for quantitative error is the truncation 
of further than nearest neighbor interactions that were 
generated by the projection to the lowest Bloch band. 
Second, the mean field approximation neglects quantum 
fluctuations of the effective spin-^ degrees of freedom. 
Inclusion of these effects certainly lead to quantitative 
changes, such as a shift of the upper critical flUing Sv^f 
at which the CDW order vanishes. Nevertheless the main 
features, including the existence of a finite window of co- 
existing CDW and superfluid orders, as well as the order- 
ing pattern of the superfluid, are expected to be robust. 

We have quantified the effect of quantum fiuctuations, 
which were neglected in the mean-field approximation by 
including spin-wave corrections. This is done using the 
Holstein-Primakoff approach (see for example Ref. 33), 
where we expand the Hamiltonian to quadratic order in 
the Holstein-Primakoff bosons. The reduced value of the 
superfiuid order parameter due to this effect is shown as 
a star in Fig. 8 for one point in the phase diagram (at 
<5j/ = 0.15). 



E. Exact diagonalization 

To support the above predictions, in particular the ex- 
istence of a supersolid phase at densities v > 1/9, we 
carry out an exact diagonalization (ED) study of the orig- 
inal kagome lattice model (1) on finite clusters. 

What are we looking for? We know from exact con- 
siderations, that at very low density v < 1/9 the sys- 
tem is infinitely compressible since filling in nonoverlap- 
ping localized flat band states does not incur an interac- 
tion energy cost. The effective theory and in particular 
the mean-fleld results imply that above the magic filling 
i^ > 1/9 one obtains a compressible state. Doping extra 
particles into interstitial sites of the CDW state incurs a 
non- vanishing interaction energy cost. This prediction is 
not obvious and is worth verifying with an exact calcu- 
lation. 



For instance, if we were to fill the lattice with fermions 
we could go up to to 1/3 filling (a full flat band) 
of polarized fermions before a non vanishing inverse- 
compressibility would set in. It is therefore natural to 
ask if an exotic, fermionized ground state of bosons could 
mimic the fermion band flUing and thus be favorable to 
the compressible supersolid state. This can be checked 
via the calculation of the compressibility as a function of 
particle doping Sv. 

Another question raised by the analytical approach 
above is about the fate of the CDW ordering in the doped 
state. From the calculation of the density-density corre- 
lation function on small clusters we can see if the corre- 
lations at the CDW wave vector Gi = ^(1/3, -l/\/3) 
are indeed enhanced as compared to a uniform state. 

We use the ALPS Lanczos application'^^ on clusters 
with a size of 27 sites corresponding to 3 x 3 unit cells of 
the kagome lattice. This corresponds to maximal dimen- 
sion d = 807885 of the Hilbert space. We perform the cal- 
culation with periodic boundary conditions shifted by a 
phase ±7r/2 along both x and y directions, which is equiv- 
alent to inserting a quarter of a flux quantum through 
each of the two holes of the torus. The twist serves to 
eliminate spurious zero energy states which wind around 
the torus. ^^ The speciflc gauge choice used to implement 
these boundary conditions is illustrated in Fig. 9 and ex- 
plained in the caption. 

The CDW order on the flnite kagome lattice is ex- 
tracted from the density-density correlation function 
Xij = {viVj). Speciflcally, we deflne the "staggered mo- 
ment" : 
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(30) 



The normalization factor A/", ensures that the staggered 



^i7r/2+i7r/2 




FIG. 9. (color online) Boundary conditions for the finite 
kagome lattice in the presence of a magnetic flux "1>± =±7r/2 
through the openings of the torus. The hopping over a blue 
bond (dashed) is associated with a phase e™ " (the bonds are 
directed as indicated by the arrow). Particles hopping over 
green bonds (dashed-dotted) bonds acquire a phase e~"^". 
The red (dash-dot-dotted) bond connecting site 9 and 20 
winds around both circumferences and picks up a phase —1. 
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10. The compressibility as a function of the density 
Vc- The mean- field compressibility is given by the 



FIG 

Su = u 

solid line while the dots indicate the results of numerical ex- 
act diagonalization. The non- vanishing compressibility of the 
CDW state &t v = Uc {5v = 0) is due to finite size of the clus- 
ter. Compressibility a,t i/ > Vc is due to the condensate frac- 
tion. The fact that the mean-field compressibility is higher 
than the exact diagonalization result is consistent with the 
higher condensate density found in the mean-field theory. 



moment is 1 in the perfect close-packed CDW state il- 
lustrated in Fig. 1. 

Fig. 8 presents the calculated CDW staggered moment 
Xsta.s ^s a function of the doping 5v away from the magic 
filling Vc — 1/9. We find Xstag = 1 at j/ = t^c as it should 
be and the moment falls off at higher filling similarly to 
the mean-field result. Quantitatively however, the CDW 
moment from ED is significantly larger than the mean- 
field moment. This may seem unusual if we compare for 
example to the Heisenberg antiferromagnet, where quan- 
tum fluctuations beyond the mean-field theory can only 
reduce the staggered moment. In our case there is the 
competing superfluid order parameter that facilitates the 
depletion of the staggered moment even at the mean-field 
level, however. Quantum fluctuations are generally more 
effective at suppressing the continuous superfluid order 
parameter then the discrete CDW. Therefore by deplet- 
ing the competing superfluid order, quantum fluctuations 
indirectly act to enhance the CDW order. Indeed we 
see in Fig. 8, that the superfluid order, measured by the 
square root of the condensate occupation number, is sig- 
niflcantly suppressed in the ED calculation as compared 
to the mean-field result. Note, that we projected the 
ED results into the flat band to make connection to the 
mean-fleld results. However, the effect of this projection 
is negligible, as we explicitly show below. 

When there are two or more competing order param- 
eters looms the danger of phase separation. To rule out 
this scenario we shall compute the compressibility n from 
the ED results using 

1/k{N) - Ea{N + 1) - 2Eo{N) + Eo[N ~ 1), (31) 

where Eq{N) is the ground state energy of N particles. 
An instability to phase separation would be signaled by 
appearance of negative compressibility. From the physics 
of the flat band below the critical doping Vc — 1/9 we ex- 
pect the system to be infinitely compressible, or k^^ = 0, 
in this regime. The CDW state itself at i^ = 1/9 is in- 
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U = Q.lt 


t/ = 0.3i 


C/ = 0.5i 


1 
2 


0.2 % 
0.4% 


1.6% 
3.0% 


3.7% 
6.7% 



TABLE I. Amplitudes of single particle states in higher bands 
extracted from the kinetic energy in the exact diagonalization 
study (see text). The low occupation of the higher bands 
justifies the effective Hamiltonians derived by projection on 
the fiat lowest band. 



compressible due to the energy gap of order U to ad- 
dition of one extra-particle on an interstitial site. The 
danger of phase separation lies when we add a finite den- 
sity of extra particles to the CDW. However the results 
shown in Fig. 10 imply a state with positive and finite 
compressibility. Moreover we see a similar trend in both 
the mean-field and ED results, both showing an increase 
in compressibility with increasing density away from the 
critical density v^ = 1/9. This gives further support to 
the mean-field picture over both phase separation and the 
more exotic possibility in which the bosons are fermion- 
ized and thereby gain high or infinite compressibility up 
to filling V — 1/3. The non-zero compressibility aX 5v — Q 
is a finite size effect due to the discreteness in the change 
of density as we increase particle number on the small 
lattice. 

Finally, we use the ED results to asses the validity of 
the projection to the lowest band. We calculate the con- 
tribution of the kinetic energy to the chemical potential 



/iM„ = ^.„(Af + l)-£;M„(^), 



(32) 



where £'ki„(A^) denotes the expectation value of the ki- 
netic energy in the ground state with N particles. To 
quantify the overlap with higher bands we take a simple 
trial state for a single quasi-particle 1-0) = cos(??)|0) -I- 
sin('i9)|l) where |0/1) denote arbitrary states in the lowest 
and the first excited band, respectively. With a finite-size 
band gap A « 0.2 i due to the twisted boundary condi- 
tions we can relate the amplitude in the first band to the 
measured kinetic chemical potential 



sin2(i9) 



A*i. 



„-2i 



(33) 



The results are given in Tab. I and we see that for all 
values of U and 5N the weight is in the range of a few 
percent. 



IV. EXPERIMENTAL REALIZATIONS 

We shall consider two types of systems that can exhibit 
the physics of bosons in fiat bands discussed above. The 
first class of systems, briefly mentioned in the introduc- 
tion, are ultracold bosons in optical lattices. The other 
class includes frustrated antiferromagnets in a large mag- 
netic field. Below we discuss the considerations relevant 
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for the realization of the flat bands in these systems and 
to the experimental detection of the emergent phases. 



A. Cold atoms 

Optical lattice potentials, formed by standing waves 
of laser light, provide the natural tool to realize inter- 
esting Hamiltonians of interacting lattice bosons. ■^^'■^''' In 
particular, there have been several suggestions on how 
to produce a lattice potential with kagome geometry. 
A direct implementation, proposed by Damski et al}'^ 
requires nine pairs of counter-propagating laser beams. 
Making use of spin dependent lattice potentials it is pos- 
sible to reduce the number of counter-propagating beams 
to three. ^^ This method also allows to construct an ef- 
fective sawtooth geometry.''*' Very recently optical lattice 
potentials were produced by an entirely different scheme 
by Bakr et al.^'' In the new scheme the optical potential 
is projected as a hologram on the ultracold atomic gas, 
allowing to easily produce any desired lattice geometry. 

Another challenging aspect of the experimental imple- 
mentation is how to obtain frustrated hopping of the 
bosons on the lattice. Usual quantum mechanical tun- 
neling of bosons between the wells naturally gives rise to 
a negative hopping matrix element, which does not lead 
to frustration. However, several clever schemes have been 
proposed to invert the hopping matrix element and thus 
produce the desired frustration. One scheme proposed 
by Eckardt et. al.^^ achieves this by introducing a fast 
time-periodic lattice acceleration. With proper design of 
the time dependent acceleration, the resulting Floquet 
Hamiltonian has an inverted sign of the hopping. The 
general idea of tuning of the effective hopping with dy- 
namical modulation'^^'^'^ was successfully demonstrated 
in experiment.'*^ 

A second scheme to invert the hopping is tailored to 
implementations of the lattice with spin dependent op- 
tical potentials. The phase of the tunneling matrix ele- 
ments can then be tuned by controlling the phase of a 
laser field that induces transitions between the two rele- 
vant spin states.''^ 

A third and completely different route for achieving 
the frustrated hopping relies on a mapping between the 
low-energy dynamics of a target Hamiltonian H and the 
dynamics in the highest energy manifold of —H.'^^ In our 
case H is the desired Hamiltonian, which describes repul- 
sively interacting bosons with frustrated hopping and we 
are interested in the low-energy dynamics on the lowest 
band (i.e. the flat band). While this is hard to produce, it 
is rather straightforward to realize —H, the un-frustrated 
model with attractive interactions. If this system can be 
initialized with all particles in the flat upper band, then 
the ensuing dynamics will be identical to the low-energy 
dynamics of the frustrated model H. Decay of particles 
from the flat band to the lowest band is exponentially 
suppressed in the large parameter hujQ/t where Hluq is 
the gap to higher Bloch bands. 



Given an experimental realization of the model (1) of 
frustrated bosons on the kagome or sawtooth lattice, the 
ground state and low-energy properties can be measured 
with standard probes. Here we concentrate on the time 
of flight expansion image which gives information on the 
ground state correlations. We note that the excitation 
spectrum can be studied using Bragg-spectroscopy'*'' or 
lattice modulation spectroscopy.'*'*'*^ 

The time of flight image lends direct access to the 
ground state momentum distribution (rik) — (fe^^k) and 
the noise correlations in the image reflect the momen- 
tum space correlations: 5(k, k') = (nk?ik') — (?^k)(?^k')■'*^ 
As we show below, these two observables complement 
each other to give detailed information on the emergent 
ground state of the frustrated bosons. 

We start by considering the momentum distribution in 
the pure CDW state on the kagome lattice a,t ly — 1/9. 
The delocalization of bosons inside hexagons gives rise to 
distinct structure in (rik), albeit no sharp delta- function 
features in absence of long-range superfluid order. To ob- 
tain the structure we rewrite the momentum distribution 
in a real space representation 



('^k) 






,jk-(xi 



--^^(blb, 



(34) 



where M is the number of kagome lattice points. Now 
using the facts that bj annihilates the CDW state unless 
j belongs to one of the occupied hexagons and that the 
off-diagonal correlations are non-zero only within a single 
hexagon we can rewrite (34) as 



K)cow = ^EE^*'^'^""°"'""'^(^"^"') 

i_^g»k.(r„-r„,) (_!)"+"' (35) 



54 



where Ri are the centers of the occupied hexagons, a la- 
bels the site of a hexagon and the vectors r^ give the lo- 
cations of the hexagon points relative to its center. Plug- 
ging in these vectors finally gives the simple expression 



("k)c 



27 



VSk^ 



which is plotted in Fig. 11(a). It is interesting to note 
that this function has zeros on the reciprocal lattice vec- 
tors of the kagome lattice (more precisely those of the un- 
derlying triangular Bravais lattice). The maxima of the 
function are located on the reciprocal lattice of the emer- 
gent CDW structure, but not on those that also belong 
to the reciprocal of the original lattice where (rik) = 0. 
Note that the calculated momentum distribution does 
not decay at large wave- vectors k. This is because we 
took lattice wave-functions to be strictly confined to lat- 
tice points. In reality the momentum distribution will 
have a decaying envelope, reflecting they are composed 
of local wave functions with flnite extent 
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(a) (nk) 



(b) (nk) 



(c) {nskn_sk) — {n5k){n_5k) 




FIG. 11. Predictions for time of flight experiments, including momentum distribution and noise correlations for the different 
phases of the kagome system, (a) The momentum distribution (rik) in the CDW state eit v = Vc- The smooth structure stems 
from the off-diagonal correlations inside hexagons. Destructive interference due to the resonating hexagon wave-function leads 
to vanishing of (rik) on the reciprocal lattice vectors of the original kagome lattice (marked by yellow crosses). The (yellow) 
hexagon marks the boundaries of the first Brillouin zone, (b) In the supersolid state at v > Uc, sharp delta-function peaks 
appear in (rik) in addition to the smooth background contribution of the CDW. These condensate peaks appear at the K and 
K' points and points related to them by the CDW wave-vector. However, points in the reciprocal lattice of the original kagome 
lattice remain exact zeros of the momentum distribution. In panel (c) we show the noise correlations for the CDW state as a 
function of the relative momentum (5k = k — k' for fixed zero average momentum. The delta peaks on the reciprocal vectors 
of the CDW lattice (peaks away from the yellow crosses that indicate the underlying kagome lattice) indicate the nontrivial 
ordering in an emergent lattice structure. 



To obtain the approximate momentum distribution in 
the supersolid phase at filling z^ > 1/9 we apply the pre- 
scription (34) to the mean-field wave-function 



*({i?„ VJJ) ) = n (sin(i^j/2) + e^"^^ cos(^,/2)w] 



Ref. 46) 



Here j runs over unit cells, or all hexagons in the kagome 
lattice and Wj creates the resonating state localized on 
hexagon j (20). As shown in section HID, -dj and Lpj 
form a three-unit cell periodic structure which amounts 
to a superfluid order parameter at the K and K' points of 
the original Brillouin zone in addition to the CDW order. 
Because of the off-diagonal long-range order, delta func- 
tion peaks appear in (rik) on top of the smooth structure 
caused by the hexagon CDW. Figure 11(b) shows the 
condensate peaks which appear at the K, K' points of 
the kagome Brillouin zone and at other points connected 
to them by the wave-vectors of the CDW. However, the 
form factor controlling the relative strength of the dif- 
ferent peaks is non trivial and most interestingly implies 
that some peaks have zero weight. Specifically, the recip- 
rocal lattice vectors of the kagome lattice remain zeros of 
(rik) although they are also connected to the K and K' 
points by CDW wave-vectors. The zeros are the result of 
destructive interference from the staggered phases of the 
resonating hexagon states. 

The noise correlations in the CDW state are computed, 
similarly to the momentum distribution, by expressing 
them in terms of the real space lattice operators (see 



^(k k') — y^ g*k-(x,-Xj)+ik'-(xi/-Xj-,) 
i'j' 

(5,,{b\b,,) + {h\by){hlh,)). (36) 



X 



The fact that {b\bj) is non zero, in the pure CDW state, 
only if i and j belong to the same occupied hexagon, 
allows us to write t/(k, k') in terms of the sites of a single 
hexagon 

g(k, k') - <55k,Q/(K, ^k) + 5,-k,G/'(K, 5k) 

/(K,Jk)=- ^ (_;L)a+a g»K.(r„-r„,)+«'5k .^ 



where 5k = k - k' and K = (k -f k')/2, Q denote the 
reciprocal lattice vectors of the original kagome lattice, 
and G are the wave-vectors of the CDW order. 

Note the sharp peaks in the noise correlations, which 
appear when k — k' is equal to a reciprocal lattice vector 
of the emergent CDW structure. This is different from 
a Mott insulating state, where lattice symmetry is not 
broken and therefore peaks are observed when k — k' is 
equal to a reciprocal lattice vector of the optical lattice 
potential."'^ The noise correlations thus provide a clear 
signature of the CDW order. 

The form factor which controls the relative strength of 
the peaks at different reciprocal lattice vectors, depends 
on both the relative momentum 5k and the average mo- 
mentum K. The dependence on the average momentum 
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is determined by the inter-site coherence within a unit 
cell of the CDW, information that is already available in 
(?ik) • It is therefore useful to fix K, as done in Fig. 11(c), 
in order to isolate the dependence of the form factor on 
(5k and thereby obtain information on the density profile 
within a unit cell. 



B. Quantum magnets 

Due to the analogy between bosons and quantum spins, 
the question of Bose condensation in flat bands has a 
lot in common with the problem of frustrated quantum 
magnets. ^^^ The frustrated geometry gives rise to an ex- 
tensive degeneracy that inhibits formation of broken spin 
symmetry, just as the flat bands inhibit Bose condensa- 
tion. 

In fact, there is a direct mapping between antiferro- 
magnets, which are nearly fully polarized by an external 
magnetic field, and the model (1) of weakly interacting 
lattice bosons. A spin flip in the otherwise polarized state 
behaves as a bosonic excitation, known as a magnon.'*^''*^ 
The magnon hops between nearest neighbor sites with a 
positive matrix element due to the antiferromagnetic ex- 
change, whereas single ion anisotropy terms oc {Sf)"^ in 
magnets with spin S >1 translate to on site interactions 
between magnons. 

Decreasing the magnetic field is like tuning the chem- 
ical potential of the magnons and can lead to Bose 
condensation. ^'^"^^ At very high magnetic fields compared 
to the Curie Weiss temperature of the antifcrromagnet, 
the spins are fully polarized and magnons are gapped. 
However, below a critical field the gap closes and the 
magnons are expected to condense giving rise to long- 
range spin correlations in the spin axes perpendicular to 
the magnetic field. Such a transition has been observed 
recently in the spin ladder system (Hpip)2CuBr4.^^ In 
the spin ladder the magnon band has a unique min- 
imum energy. The approach developed in this paper 
becomes useful when the lattice geometry is sufficiently 
frustrated to induce a flat magnon dispersion, ^^ which 
inhibits straightforward magnon condensation. 

An interesting system from this viewpoint is the 
organic material m-MPYNN-BF4,-'^^'-'^^ which forms a 
kagome structure of effective spin-1 degrees of freedom. 
Fortunately it also has a low Curie- Weiss temperature 
(©cw ~ 3K), which allows to completely polarize the 
spins with reasonable magnetic fields. Indeed, NMR spin 
relaxation measurements indicate the presence of a spin 
gap above a field /ic2 ^ 4T.^'* The close packed hexagon 
lattice should show up as a magnetization plateau at 8/9 
of full polarization at magnetic fields slightly below the 
saturation field. Our results pertain to the nature of the 
ground state with magnetization below 8/9 per spin. 

The detailed ordering structure can in principle be ob- 
tained from elastic neutron scattering experiments which 
measure the spin correlation function. In particular, the 
X — y component of the static structure factor {S^ S~) 



can be mapped directly to the momentum distribution of 
the bosons discussed above and shown in Fig. 11(b) (see 
App. B for the details of the spin-boson mapping). The 
z-component of the spin structure factor can be directly 
mapped to the Fourier transform of the density correla- 
tions. It will show peaks at the CDW wave-vectors mod- 
ulated by a form factor which depends on the density 
profile within a unit cell. In that sense neutron scatter- 
ing will give information similar to the noise correlations 
in time of light images of ultracold atoms at a fixed value 
of the center of mass momentum (See Fig. 11). 

Finally, we note that in order to apply the results of 
our study to m-MPYNN-BF4, it may be necessary to 
take into account the structural distortion that occurs 
below Ts « 130 K. This distortion leads to anisotropy of 
the magnetic exchange along different directions of the 
kagome lattice. ^^ As long as the anisotropy is weak com- 
pared to the effective magnon interactions then the use 
of the effective Hamiltonian restricted to the lowest band 
is still justified, albeit with spatially anisotropic interac- 
tions. 



V. CONCLUSIONS AND OUTLOOK 

We have developed a general scheme to derive low- 
energy effective Hamiltonians for bosons on the ffat Bloch 
bands of frustrated lattices. The essential idea is to 
project the Hamiltonian on the flat band by formulat- 
ing the interaction term using a complete basis of lo- 
calized Wannier states of the flat band. The resulting 
model contains both diagonal interaction terms and off- 
diagonal (assisted) hopping terms, all proportional to the 
interaction parameter in the original model. Most im- 
portantly, the effective model is in most cases no longer 
highly frustrated and therefore amenable to analysis by 
mean-field theory and other standard theoretical meth- 
ods. Such analysis is directly useful in the description of 
ultracold atoms in optical lattices with highly frustrated 
hopping and of frustrated quantum magnets in high mag- 
netic fields, where magnons behave as dilute bosons. 

Using the projection scheme, we investigated the Bose 
Hubbard model on the frustrated sawtooth and kagome 
lattices. In both cases it is possible to construct exact 
many-body ground states composed of nonoverlapping 
configurations of localized states, which lie entirely in 
the ffat band^"'^^ (see Fig. 2). The existence of a large 
number of such configurations prevents any ordering at 
zero temperature. Above the critical lattice filling Vc^ at 
which the localized states form a closely packed density 
wave, the exact construction fails and our approach is 
useful in identifying the emergent ground state and low- 
energy excitations. 

In the one-dimensionalsawtooth lattice we show that 
the CDW is immediately destroyed by proliferation of 
domain walls upon increasing the density from the magic 
filling Vc — 1/4. Universal properties, e.g., the power- 
law peaks in the static structure factor, at filling v = 
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1/4: + Sv are described by the long- wavelength theory of 
the coniinensurate-incommensurate transition. 

In the kagome lattice, the CDW remains stable upon 
increasing the density from Vc = 1/9. The added bosons 
hop between interstitial sites of the CDW and form a 
condensate at the K and K' points of the Brillouin zone 
of the kagome lattice. We considered the signatures of 
this phase in a realization with ultracold atoms in optical 
lattices. The predicted time of flight image is shown in 
Fig. 11(b) shows the sharp peaks indicating the presence 
of a condensate at the right wave- vector. Most strikingly 
however, as a result of destructive interference from the 
localized hexagon wave-functions, we find zeros of the 
momentum-distribution on the reciprocal lattice vectors 
of the original kagome lattice. 

A general feature of Bose condensates in flat bands of 
frustrated lattices is that condensation is entirely inter- 
action driven. It is a result of the emergent dynamics 
induced inside the flat band by the weak interactions, 
which is directly captured in the projection scheme. In 
particular both the condensation energy and superfluid 
stiffness are proportional to the Hubbard interaction U 
and not to the kinetic energy as in usual condensates. 
In a similar way, due to effective interactions generated 
in the low-energy model, charge density wave states and 
supersolids are stabilized even in absence of long-range 
interactions in the microscopic model. 

It would be interesting to apply the ideas developed in 
this paper to frustrated antiferromagnets such as the or- 
ganic salt m-MPYNN-BF4. Such analysis should include 
coupling to the anisotropic lattice distortion. Further- 
more, nearest neighbor magnon interactions, which are 
naturally generated by the antiferromagnetic magnetic 
exchange, are expected to give rise to more magnetiza- 
tion plateaux. An intriguing question then is how the 
supersolid phase we find connects with the plateaux at 
lower magnetization? Using the projection method to 
investigate the compressible phases, which connect be- 
tween lower plateaux, and perhaps finally with the spin 
gapped phase at zero field, offers a new viewpoint on the 
problem of frustrated quantum magnets. 
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FIG. 12. (color online) Wannier functions (probability distri- 
bution) as a function of position on the sawtooth lattice with 
t' /t = \pl. The excited band has an oscillatory slowly decay- 
ing tail, while the flat-band Wannier functions decay with a 
localization length ^ « log(2.15)a. The inset shows that lo- 
cally, the Wannier function of the flat band is almost identical 
to the strictly localized V states indicated by the fllled circles. 



Appendix A: Details of the projection to the flat 
band 



The explicit expressions for the Wannier functions of 
the flat band of the sawtooth lattice can be written in 
terms of complete elliptic integrals of the first and second 
kind, K(m) and E(m), respectively. 

Specifically, the amplitudes used in the derivation of 
the effective Hamiltonian (10) are given by 



wb{0) = 
wsia) = 

wa{0) = 



30? 
2 

2 1 



K{2/3), 

{3E{2/3) + 2K{2/3)), 



Stt 



{3E{2/3)-K{2/3)), 



(Al) 
(A2) 

(A3) 



2 1 



WA{a)^^--{3K{2/3)-5Ei2/3)). (A4) 

The resulting interaction constants are of the form: 

F/U w 8wl{0)w%{a) + 8w\{0)w\{a) + 4w\{0) 

w 0.112, (A5) 

/""/C/ « 4w|(a)ws(0) + 4w\{0)wAia) + Awl{a)wAiO) 

w -0.025, (A6) 

4wl{a)wl{0)+8w\iO)wAia) 

-0.011, (A7) 

w%{0) + 2wB{af + 2w\{0) + 2w\ia) 

0.40. (A8) 



r^/u 
ii/u 



Fig. 12 shows the spatial dependence of the Wannier 
functions of both the lower flat band and the upper band. 
Note the strong localization of the Wannier function of 
the flat band. A zoom-in on the short range structure 
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FIG. 13. (color online) Wannier function for the flat band in 
the kagome lattice, (a) Wave-function amplitudes shown on 
the two-dimensionallattice, where the size of the dots repre- 
sents the weight and the color (shading) the denotes the sign 
+ or — of the wave-function. The definitions of the weights 
wi . . .W5 are indicated, (b) Algebraic decay of the Wannier 
function as ~ l/l?"! along the high symmetry direction marked 
by the dashed line in (a). 



of this function on-site and on nearest neighbors shows 
that it is shows that it is almost identical to the strictly 
localized V states of the sawtooth lattice. 

Fig. 13 shows the spatial behavior of the flat band 
Wannier functions of the kagome lattice as defined in 
(22). Due to the touching band, they do not fall of ex- 
ponentially but decay algebraically, cf. Fig. 13(b). ^^ 

Coefficients in Hamiltonian (23) result from the over- 
lap of Wannier functions at different sites. The explicit 
expressions for the first few terms are given by 



U/U^ 


^ Qwl « 0.144, 




(A9) 


F-/U^ 


a Wq + Awlwl + 2wlwl ? 


a 0.028, 


(AlO) 


il^iu. 


a 2wlwx « 0.009, I^y -- 


= /r/2, 


(All) 


i^y/u. 


a 2wlw2 ~ 0.003, ll^ -- 


- ^3 V2 


(A12) 



which includes anisotropic magnetic exchange interac- 
tions as well as a single ion anisotropy D. The latter 
is relevant for S > 1/2 where two magnons can reside on 
the same site; we concentrate on this regime in the fol- 
lowing. This Hamiltonian has been studied extensively 
at small magnetic field h and S — l.^^^^o 

Here, we focus on the large field limit. The fact that 
spins are nearly polarized in the z direction suggests to 
construct a bosonic representation of the small fluctua- 
tions (magnons) about the fully polarized state using the 
Holstein-Primakoff expansion. '^'^ To quadratic order in 
the magnon occupation the mapping between the spins 
and bosons is: 



^.^-^fl-%)&. 



b]b, 



si^bi^2syi-^ 



(B2) 



{S-blb: 



Expressing Hamiltonian (Bl) using the Holstein- 
Primakoff bosons and keeping terms to leading order in 
1/ S we finally obtain the bosonic Hamiltonian 



H^JS^[b\b^+}i.c]+Y^ Me„.z>, + D^v^, (B3) 



(«j) 






bj + H.c] + XviVj 



0(1/^2 



where the amplitudes Wa are defined in Fig. 13. 



Appendix B: Spin— boson mapping 

In this appendix we describe details of the mapping 
from a spin-S' system to a bosonic Hamiltonian. We con- 
sider the following general spin Hamiltonian 






(Bl) 



Here, Oi — b\b^ denotes the density operator. The effec- 
tive chemical potential for the magnons here is given by 
;x,jj = SD + 2z\JS + h. The Bose-Hubbard model (1), 
which was considered in this paper, is only the first line of 
(B3). The second line includes additional nearest neigh- 
bor interactions and assisted hopping terms which stem 
from the antiferromagnetic exchange interactions. Note 
however, that both terms in the second line are vanish- 
ing in the exact CDW state and can be projected into 
the flat band along the same lines as the simple on-site 
interaction ex D. 
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